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Introduction 

The NASA self calibrating multiwavelength pyrometer* 1 ’ is a recent addition to the list of 
pyrometers used in remote temperature measurement in research and development. The older one- 
color, two-color and the disappearing filament pyrometers* 2 ' as well as the multicolor and early 
multiwavelength pyrometers* 3-5 ' all do not operate successfully in situations in which strong ambient 
radiation coexist with radiation originating from the measured surface. Figure 1 depicts the situation 
in question. Radiation departing from the target surface arrives at the pyrometer together with 
radiation coming from another source either directly or through reflection. U nlik e the other 
pyrometers, the self calibrating multiwavelength pyrometer can still calibrate itself and measure the 
temperatures in this adverse environment. 


Theory 


In reference 1, it was shown that voltage spectra measured by a pyrometer contain enough information 
to determine all the quantities necessary for pyrometry. The voltage spectra V(X,T(t)) are of 
wavelengths X v X,, ..X s .. X*. N>2, at times t v t 2 , .. t,, .. t n , though measured at unknown temperatures 
T(t x )*T(t 2 )* .. *T(t;>* ..*T(t n ), with unknown target surface spectral emissivity £*, unknown optical 
medium transmissivity \ between the measured surface and the pyrometer , and unknown instrument 
calibration constant g x , all these unknowns as well as the unknown temperatures can be determined. 
The bases of the self calibrating pyrometer is contained in three equations* 1 ’ (Eqns. 1-3) connecting 
two 2 wavelengths, labeled Xr and X T . 
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With the determination of a single number, g R e R T R , all the other g^t, and temperatures T(t) are 
automatically determined. g R e R x R is determined using a least squares variational method, beginning 
by assuming a value for g R e R x R , plotting the quantity on the left hand side of Eqn. 1 vs the quantity 
inside the curly bracket on the right hand side to obtain a straight line, the slope gives Do this 

for all I*R to determine g^ft, and use them to calculate the average temperatures T(t). 


In the application of this multiwavelength pyrometer to the geometry of figure 1, because of the 
presence of the ambient radiation, the voltage spectrum, now denoted by S(X,t) is given in Eqn. 4. 
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I(X,t) is the spectral intensity of the ambient radiation and I’(X,t) is a transformation obtained from 
it. If I(X,t) and hence I’(X,t) is slowly and randomly varying with time, then the quantity V(X,t), in 
Eqns. 1 to 3, intrinsically related to the temperature of the measured surface can be recovered and 
used to determine the relevant g, e, x and hence the temperatures T(t) by the following steps: 

SI: Beginning with t=0, a difference spectrum DiX^.O) is formed according to Eqn. 5, between the 

first and second measurement. 

D(X, t lt t 0 =0) =(S(X, t x ) ~S(X, t 0 =0) ) 
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More generally the difference spectrum is defined iteratively in Eqn. 6. 

D(X, t ± , t d . d ) = (SU, t d ) ~S(X, t ± . d ) ) 
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S2: Form the sum of the difference spectra from t=0 to t=tj according to Eqn. 7. 


£d(X, t ± , t i . 1 ) = (S(X, t d ) -S(X, t=0) ) 
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S3: Choose a reference time t,. At this time the temperature is T r . Analogous to Eqn. 7, define 

the sum of the difference spectra from t=0 to t=ty in Eqn. 8. 
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S4: The difference between Eqn. 7 and Eqn. 8 (Eqn. 9) is the sum of all the difference spectra 

between t=t r and t=tj 

i r 
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S5: This is rearranged in Eqn. 10. 

g x e x x x L(X, t d ) =g x e x x x L(X , t z ) +AS(X, t ± ) +Al'(X, t s ) (10) 

where L(X,t) is the Planck formula, and AS(X,tj), AT(X,tj) are defined as 

AS" (A , fcj) = S(X, t z ) - S(A., t z ) (11) 

AJ'U, =J'U, t ± ) -l'(X, t z ) 

S6: For constant or slowly and randomly varying I(X,t), AT(X,t)=0, Eqn. 10 becomes 

g x t x x x L(X , tj) =g x t x z x L{X , t z ) +AS(X r t ±) (12) 

g x e x 't x L(A,,t i ) is the voltage spectrum that the pyrometer would have detected at t, in the absence 
of the ambient radiation 1(1,^), and it equals g^x^IXXjt,), the voltage spectrum that the 
pyrometer would have detected at time in the absence of the ambient radiation plus AS(l,tj), 
the increase in the sum of difference spectra when the target temperature has changed from 
T(t,) to T(tj) between time t, and time t ; . The expression for V(X,t) is now given by Eqn. 13. 

V(X, fcj) =g x t x x x L(X , t z ) +A S(X, t ± ) (13) 


Experiment and Results 


The voltages V(X,t) in Eqns. 1 to 3 in the derivation of the self-calibrating multiwavelength pyrometer 
in reference 1 are now given by Eqn. 13. By direct substitution and minor modification, using 
difference spectra obtained in steps (SI to S6) are processed to determine the unknown temperatures 
and all the necessary calibration constants. Substituting Eqn. 13 in Eqn. 1 and rearranging gives* 
Eqn. 14. 
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Referring to Eqn. 14, assign a value for g R e R x R and a value for T r , evaluate the quantity inside the 
curly bracket on the right hand side and use it to plot against the quantity on the left hand side of 
it. A straight line is obtained. Standard least squares method is used to determine its slope, which 
is g[£jT[. Another least squares procedure is used to determine g R £ R x R and T r as follows: 


1) 

2 ) 

3) 

4) 

5) 

6 ) 


7) 


Choose a value for T r . 

Choose a wavelength for and a wavelength for X, and a value for g R £ R t R . 

Plot the data according to Eqn 14 to determine the g^,!, from the slopes. 

Do so for all Ml. 

Use the so determined g^Xj to calculate the temperatures T(t) according to Eqns. 2, 3 and 13. 
Transform the spectra V(X,,t) into a single large data set (x,y). The wavelength X, is 
transformed into a/XjTCt), the transformed wavelength, the voltage V(Xj,t) is transformed by 
first dividing by g^x, and then by T(t) 5 , the 5th power of the spectrum temperature, according 
to the prescription 


i _ c i x 5 

XT(t) A T(t) s d e x -l 


(15) 


The transformed (x,y) data obey the generalized non-dimensional Planck function. The (x,y) 
data are fitted to the Planck function in Eqn. 15 by calculating the residual I, defined as the 
sum of the squares of the difference between the transformed y] and the calculated y evaluated 
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by substituting the transformed x, in the y equation in Eqn 15 for all the data. 

8) A new value for g R e R T R is selected, and steps (2 to 7) repeated. 

9) The value of g^^ that produced the least X is recorded. 

10) Choose a new value for T r , repeat steps (2 to 9). 

11) The combination of g R e R x R and T r that produced the least over all X is the correct one we are 
after. 

When g R E R x R) T r are determined, the other T(t) are also determined. In one stroke, everything that is 
needed in pyrometry for temperature measurement is determined. Temperatures at any time in the 
past or in the future are determined from each X according to Eqn 3 or by least squares curve fitting. 
As is evident, the pyrometer requires no prior calibration. The procedures (1 to 11) can be repeated 
as often as necessary during an experiment to update the self calibrating process. 

A simulation experiment of the condition in figure 1 is shown in figure 2. It consisted of a black body 
furnace and the spectrometer of a multiwavelength pyrometer; between them were placed a quartz 
lamp and a slab of transmitting material (zinc selenide, ZnSe). This transmitting material was 
polished so that it possessed optical quality parallel surfaces. It was then positioned between the 
quartz lamp and the spectrometer so that the black body radiation transmitted through it, and the 
quartz lamp radiation, specularly reflected, would arrive and mix at the pyrometer detector. .The 
ZnSe-transmitted black body radiation simulated the target surface radiation; the ZnSe transmissivity 
is unknown and simulated the target surface emissivity; the more intense, specularly reflected quartz 
lamp radiation in this simulation replaced the weaker, diffusely reflected ambient radiation normally 
would be encounter in an actual experiment. The quartz lamp emitted (1) short wavelength (X < 2.5 
pm) radiation transmitted from the filament through its envelope and (2) longer wavelength radiation 
emitted by the envelope after its temperature had increased. 

In this experiment, the spectrometer of a multiwavelength pyrometer recorded a series of spectra. A 
stable power supply generated a constant current through the quartz lamp filament to produce a 
constant radiation flux. The temperature of the black body furnace was gradually increased. The 
Klarlr body furnace temperature was measured using a type G thermocouple. The TC measured black 
body furnace temperatures are 763, 871, 984, 1101, 1224 and 1351 K. Six of the changing spectra are 
shown in figure 3. These spectra spanned the spectral region from 1.3 to 14.5 pm. It recorded the 
direct voltage output of the indium antimonide and mercury cadmium telluride detectors. 

Results The spectrum minim a in figure 3 are due to atmospheric C0 2 and HjO absorptions iq 

the optical path between the detector and the black body source. One of the spectra in figure 3 is 
chosen as the reference, it is subtracted from the other spectra, including itself. The difference 
spectra so generated are shown in figure 4. Following the analysis above, Xr is chosen to be 5 pm. 
An arbitrary initial value for g R e R x R is used, plots according to Eqn. 14 are made. Figure 5 shows the 
case for 5^=5 pm, and \=2 pm, g R e R x R = 1.6468, and T r =1224 K. It is indeed a straight line. The slope 
(i.e. gj£|Xj) of plots like this at other wavelengths are obtained using least squares method. They are 
plotted in fig. 6. These g^x, are now used to calculate the temperature of each spectrum according 
to Eqn. 2. The intrinsic voltages are given by Eqn. 13, and the calculated temperatures are shown in 
fig. 7. They are alm ost independent of wavelength. The explosion at the shortest wavelength and low 
temperatures are due to poor signal to noise there. At some wavelengths, there is no solution to Eqn. 
2 which involves the logarithm. When data at these wavelengths are excluded, averages are obtained 
according to Eqn. 3. These averages, and everything that depends on g^ R x R change when g R e R x R takes 
on different values. Of the many possible values that g R e R x R can assume, we identify the correct one 
using the least squares curve fitting procedure described in steps 1 to 11. Least squares fitting of the 
transformed data to the non-dimensional Planck formula is shown in figure 8. The agreement is 
excellent. 

Once the value of g R e R x R and T r are fixed, the values of all other g,£[X, and the temperatures of the 
other spectra are fixed as a consequence. These temperatures are 750, 864, 981, 1101, 1224 and 1354 
K. The pyrometry measured temperatures and the TC measured temperatures are compared in fig. 9. 
They differ by less then 2%. Obviously the pyrometer measured these temperatures without having 
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been previously calibrated. 

Four of the 445 wavelengths were selected for a similar analysis, a smaller (x,y) data set is now 
generated. These wavelengths were 2.4, 3, 4 and 5 pm. The value of g R e R x R that produced the least 
squares is 1.647, almost exactly the same as before. The resulting fit of the data to the generalized 
Planck function is performed and shown in figure 10. The temperatures determined were also almost 
exactly the same. The choice of these 4 wavelengths occur in the regions where the signal is strongest 
in the temperature range that measurements were made. 

Conclusion 

The multiwavelength pyrometer successfully measured the temperatures of a black body furnace 
viewed through a transparent window of unknown transmissivity, in the presence of constant or slowly 
and randomly varying strong ambient radiation. The measurement error is less than 2%. This is 
significant because the pyrometer is not previously calibrated with a standard such as a black body 
furnace. The use of only 4 wavelengths in the most optimal spectral regions produced almost identical 
results, reducing the amount of data and data acquisition time by a factor of over 100. 
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Surface 




Arrangement for Simulation Experiment. 





Figure 3 

Spectra of Simulated Quartz-heated Surface. 



Figure 4 

Difference spectra obtained by subtracting a chosen reference from the other 
spectra. 
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Figure 7 

Calculated temperatures use the determined calibration constants. 


Transformed Spectra Fitted 
to Generalized Planck Function 



Figure 8 

Least squares curve fitting of the transformed data to the non-dimensional Planck 
function. 
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Figure 9 

Plot of pyrometer measured temperature vs thermocouple measured temperature. 


Transformed Spectra Fitted 
to Generalized Planck Function 
(Use only 4 Wavelengths) 



Figure 10 

Fit of data from 4 wavelengths to the non-dimensional Planck Equation. 
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